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In a previous Letter [Phys. Rev. Lett. 77, 4158 (1996)], 
a new correlation measure was introduced that sensitively 
probes phase space localization properties of eigenstates. It 
is based on a system's response to varying an external pa- 
rameter. The measure correlates level velocities with overlap 
intensities between the eigenstates and some localized state of 
interest. Random matrix theory predicts the absence of such 
correlations in chaotic systems whereas in the stadium bil- 
liard, a paradigm of chaos, strong correlations were observed. 
Here, we develop further the theoretical basis of that work, 
extend the stadium results to the full phase space, study the 
^-dependence, and demonstrate the agreement between this 
measure and a semiclassical theory based on homoclinic or- 
bits. 

PACS numbers: 03.65. Sq, 05.45. Mt 



I. INTRODUCTION 

The two general motivations for our investigation 
are understanding better the nature of eigenstates of 
bounded quantum systems possessing 'simple' classical 
analogs, and exploring new features of such systems' be- 
havior as a system parameter is smoothly varied. Sim- 
ple in this context refers to few degrees of freedom and 
a compact Hamiltonian. Nevertheless, the classical dy- 
namics may display a rich variety of features from regu- 
lar to strongly chaotic motion. We focus on the strongly 
chaotic limit for which semiclassical quantization of in- 
dividual chaotic eigenstates does not hold, and the cor- 
respondence principle is less well understood [fjj. Even 
though there has been some recent progress Q , it turns 
out that with a detailed understanding of chaotic sys- 
tems a statistical theory provides a well-developed, al- 
ternative approach to these difficulties. Twenty years 
ago, Berry ||] conjectured and Voros || discussed that 
in this case as h — > the eigenstates should respect the 
ergodic hypothesis in phase space, S(E — H(p, q)), as 
it applies to wavefunctions. In essence, the eigenmodes 
should appear as Gaussian random wavefunctions locally 
in configuration space with their wavevector constrained 
by the ergodic measure of the energy surface. Discussion 
of the properties of random waves and recent supporting 
numerical evidence can be found in refs. (^||. 

The second general motivation relates to a long recog- 
nized class of problems, i.e. a system's response to para- 



metric variation. Our interest here is restricted to ex- 
ternal, controllable parameters such as electro-magnetic 
fields, temperatures, applied stresses, changing bound- 
ary conditions, etc., through whose variation one can 
extract new information about a system not available 
by other means. A multitude of examples can be found 
in the literature 0. A recent concern has been univer- 
salities in the response of chaotic or disordered systems 
and statistical approaches to measuring the response [0 . 
Universal parametric correlations have been derived via 
field theoretic or random matrix methods for quanti- 
ties involving level slopes (loosely termed velocities in 
this paper), level curvatures, and eigenfunction ampli- 
tudes (^[lO). In contrast, our motivation is not the uni- 
versal features per se for they cannot tell us anything 
specific about the system other than it is, in fact, chaotic 
and/or symmetry is present. Rather we are interested 
in what system specific information can be extracted in 
the case that the system's response deviates from univer- 
sal statistical laws. The specific application discussed in 
this paper shows how one can decipher phase space local- 
ization features of the eigenstates. The theory naturally 
divides into a two-step process. One must first under- 
stand any implied limiting universal response of chaotic 
systems. Next, one must develop a theory which gives 
a correct interpretation of any deviations seen from the 
universal response. The necessarily close interplay be- 
tween theory and observation required to deduce new 
information forms part of the attractiveness of investi- 
gating parametric response. 

Taking up the first step of understanding universal re- 
sponse, an expected but rarely discussed property is the 
independence of eigenvalue and eigenfunction fluctuation 
measures jll| which is found in the random matrix the- 
ories anticipated to describe the statistical properties of 
quantum systems with chaotic classical analogs p^fD!. 
Coupled with Berry's conjecture mentioned above, these 
properties imply a 'democratic' response to parametric 
variation for an ergodically behaving quantum system. 
The perturbation connects one state to all other states 
locally with equal probability. The variation of any one 
eigenstate or eigenvalue over a large enough parameter 
range will be statistically equivalent to their respective 
neighboring states or levels. 

In a pioneering work on the ergodic hypothesis using 
the stadium billiard, now a paradigm of chaos studies, 
McDonald noticed larger than average intensities of the 
eigenstates in certain regions [H . In his thesis he states 
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that "a small class of modes (bouncing ball, whispering 
gallery, etc.) seem to correspond naively to a definite set 
of 'special' ray orbits." Heller initiated a theory concern- 
ing these large intensities when he modified the random 
wavefunction picture with his prediction and numerical 
observations of eigenstate scarring Jl5|. He derived a 
criterion for eigenstate intensity in excess of the ergodic 
predictions along the shorter, less unstable periodic or- 
bits. Scarring is thus one possible phase space 'localiza- 
tion' property of a chaotic eigenstate. Other possibilities 
result from time scales not related directly to the Lya- 
punov instability such as transport barriers in the form 
of broken separatrices , and cantori [ p7|p8| , or diffu- 
sive motion jljj). In the context of this paper, we take 
localization to mean some deviation from the ergodic ex- 
pectation beyond the inherent quantum fluctuations, and 
it creates the possibility of a non-democratic response to 
parametric variation. A perturbation could preferentially 
connect certain states or classes of states, thus leading 
to additional short-range avoided crossings or like level 
movements within a particular class, etc. 

Debate ensued Heller's work on eigenstate scarring, in 
part, because of the difficulty in quantitatively character- 
izing and predicting its extent in either a particular eigen- 
state or even collective groups of eigenstates. Judging 
from the earlier literature, it was easier to graph eigen- 
states in order to see the scarring by eye than define pre- 
cisely what it means or what its physical significance is. 
Furthermore, he linearized the semiclassical theory which 
was insufficient for a full description of scarring. We re- 
mark that recent work suggests the opposite, i.e. the lin- 
earized theory is sufficient assuming h is smaller than 
some system specific value which is 'small enough' [ p0[ . 
However, many of the experimental and numerical in- 
vestigations are far from this regime and the nonlinear 
dynamical contributions are essential for understanding 
most of the work being done. The theory incorporating 
nonlinear dynamical contributions [p|]2l|l was developed 
much later than Heller's introduction of scarring. It is 
based on heteroclinic orbit expansions for wave packet 
propagation and strength functions. Ahead, we make ex- 
tensive use of these forms to derive a semiclassical theory 
applicable to problems involving parametric variation. 

In a previous Letter |z| , one of us (ST) introduced a 
measure that very sensitively probes phase space localiza- 
tion for systems having continuously tunable parameters 
in their Hamiltonians. It correlates level motions un- 
der perturbation with overlap intensities between eigen- 
states and optimally localized wave packet states. The 
basic idea is that the wave packet overlap intensities se- 
lect eigenstates that potentially have excess support in 
the neighborhood of the phase point at the wave packet's 
position and momentum centroids. The perturbation will 
push these levels somewhat in the same direction depend- 
ing on how it is distorting the energy surface near that 
particular phase point. If the level velocities associated 
with those states have similar enough values, then sig- 
nificant non-zero correlations will result that reveal the 



localization. The measure can be used in a forward or 
reverse direction. If phase space localization is present 
in a system of interest, then it predicts experimentally 
verifiable manifestations of that localization. Conversely, 
one can first experimentally determine the level velocity - 
overlap intensity measure in that system for the purpose 
of inferring the existence and extent of localization. 

Our purpose in this paper is to give a complete account 
of that Letter, develop further the semiclassical theory, 
and explore the full phase space and h behavior of the sta- 
dium billiard, a continuous time system. In a companion 
paper immediately following this one, we give the the- 
ory for quantized maps (discretized time) p3| . The next 
section introduces strength functions and a new class of 
correlation coefficients. Section III utilizes ergodicity and 
random wave properties to motivate the introduction of 
random matrix ensembles. The ensembles describe the 
statistical properties of chaotic systems in the h — > 
limit. The correlation measures vanish for these ensem- 
bles indicating the absence of localization and universal 
response to perturbation (i.e. parameter variation). Sec- 
tion IV gives the semiclassical theories of level velocities, 
and overlap intensity-level velocity 
We finish with a full treatment 



strength functions, 
correlation coefficients, 
of the stadium billiard and concluding remarks. 



II. PRELIMINARIES 

Consider a quantum system governed by a smoothly 
parameter-dependent Hamiltonian, H(X) with classical 
analog H(p, q;A). We suppose that the dynamics are 
chaotic for all values of the A range of interest, and sup- 
pose the absence of symmetry breaking. Then the expec- 
tation is that all statistical properties are stationary with 
respect to A. Without loss of generality, we also assume 
the phase space volume of the energy surface is constant 
as a function of A. This ensures that the eigenvalues 
do not collectively drift in some direction in energy, but 
rather wander locally. We use the same strength function 
Heller employed in his prediction of scarring |]l5| except 
slightly generalized to include parametric behavior; 

1 r°° 

Sa (E, A) = — — / dt e^He-^^la) 



= Tr[p a S(E - H(X))] 

= Y J P<*n{\)8{E - E n {\))- 



(1) 



Pan (A) = \{a\E n {\)){ 



where p a — \a){a\. S a (E, A) is the Fourier transform of 
the autocorrelation function of a special initial state \a) 
of interest. Ahead S a (E, A) will denote the smooth part 
resulting from the Fourier transform of just the extremely 
rapid initial decay due to the shortest time scale of the 
dynamics (zero- length trajectories). We will take |a) to 
be a Gaussian wave packet because of its ability to probe 
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"quantum phase space," but other choices are possible. 
Say momentum space localization were the main inter- 
est, the natural choice would be a momentum eigenstate. 
| a) can be associated with a phase space image p Q (p, q) 
of Gaussian functional form using Wigner transforms or 
related techniques. p a (p, q) turns out to be positive def- 
inite and maximally localized in phase space, i.e. it oc- 
cupies a volume of h d . 

For a fixed value of the parameter, an example strength 
function is shown in Fig. ([!]). If the wave packet is cen- 
tered somewhere on a short periodic orbit, large ampli- 
tudes necessarily indicate significant wave intensity all 
along the orbit as seen in the inset eigenstates. This 
behavior cannot, a priori, be stated to be obviously in 
violation of the quantum statistical fluctuation laws even 
if it appears so. That remains to be determined. With 
the inclusion of parametric variation, the eigenvalues of 
a chaotic system are supposed to move along smoothly 
varying curves of the type shown in the upper square of 
Fig. (§). Many of the previous studies of parametric vari- 
ation focussed on the properties of such level curves. A 
great deal is known about the distribution of level veloci- 
ties p^p5| , the decay of correlations in parametric statis- 
tics ]10| , |26[ , the distribution of level curvatures p7] , 
and the statistics of the occurrences of avoided cross- 
ings MM. 
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FIG. 1. Strength function for the stadium billiard. The 
Gaussian wave packet is centered in the stadium with mo- 
mentum directed towards the end cap. The large intensities 
are where the scarring occurs. 



A new class of statistical measures can be defined that 
cross correlate intensities with levels. The most evident 
examples are the four correlation coefficients involving 
both level curves and eigenstate amplitudes that can be 
defined from the following quantities: (i) the level veloci- 
ties, 8E n (X)/dX, (ii) level curvatures, d 2 E n (X)/dX 2 , (iii) 
overlaps, p ani and (iv) overlap changes, dp an /dX. The 
most important is the overlap intensity-level velocity cor- 
relation coefficient, C Q (A), which is defined as 



C a {X) 



dE n (X) \ 
fa™ d\ J E 



(2) 



where a 2 a and o\ are the local variances of the over- 
laps and level velocities, respectively. The brackets de- 
note a local energy average in the neighborhood of E. It 
weights most the level velocities whose associated eigen- 
states possibly share common localization characteristics 
and measures the tendency of these levels to move in a 
common direction. In this expression, the phase space 
volume remains constant so that the level velocities are 
zero-centered (otherwise the mean must be subtracted), 
and C a (X) is rescaled to a unitless quantity with unit 
variance making it a true correlation coefficient. The set 
of states included in the local energy averaging can be 
left flexible except for a few constraints. Only energies 
where S a (E, A) is roughly constant can be used or some 
intensity unfolding must be applied. Also, the energy 
range must be small so that the classical dynamics are 
essentially the same throughout the range, but it must 
also be broad enough to include several eigenstates. 




We now superpose the strength function overlap inten- 
sity information on Fig. (^|) in the lower square as vertical 
lines centered on the levels; the lengths are scaled by the 
intensities (3-D versions of this figure turned out not to 
be very helpful). By considering the full strength func- 
tion and not just the level curves (i.e. density of states), 
the eigenstate properties can be more directly probed. 



FIG. 2. Illustration of ergodic behavior. The upper square 
shows how the energy eigenvalues move as a function of A. 
The lower square is a graphical representation of S a (E,X). 
Each small line segment is centered on an eigenvalue and its 
lambda value. The heights are proportional to the overlap 
intensity with a wavepacket. The level velocities and overlap 
intensities were produced using a Gaussian orthogonal ensem- 
ble. 
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C a (A) thus has a simple form and the additional ad- 
vantage of involving quantities of direct physical inter- 
est. Level velocities (curvatures also) arise in thermody- 
namic properties of mesoscopic systems |3^ j, and over- 
lap intensities often arise in the manner used to cou- 
ple into the system |3^ j. It is the most sensitive mea- 
sure of the four possible combinations, the others being 
the intensity-curvature, intensity change - curvature and 
intensity change - level velocity correlation coefficients. 
The first two are far less sensitive measures of eigen- 
state localization effects, even though curvature distribu- 
tions are affected by localization because of the relative 
rareness of being near avoided crossings where curvatures 
are large. The last shows no effect since intensities will 
change whether the level is moving up or down. These 
three measures will not be considered further in this pa- 
per, but we did calculate them to verify their lack of 
sensitivity. 

III. ERGODICITY, RANDOM WAVES, AND 
RANDOM MATRIX THEORY 

Semiclassical expressions for wavefunctions have the 
form 

*(x) = M*) exp (%S n (x)/h - iu n ir/2) (3) 

n 

where S n (jx.) is a classical action, v n is a phase index, and 
A n (x) is a slowly varying function given by the square 
root of a classical probability. The classical trajectory 
underlying each term arrives at the point x with mo- 
mentum, p„ = VS n (x). For a chaotic system, a complete 
theory leading to an equation of the form of Eq. (^) does 
not exist [fij. Nevertheless, Berry [[| conjectured that 
for the purposes of understanding the statistical prop- 
erties of chaotic cigenfunctions, the ergodic hypothesis 
implies that the true eigenfunction will appear statisti- 
cally equivalent to a large sum of these terms each ar- 
riving with a random phase (since each wave contribu- 
tion extends over a complicated, chaotic path). For sys- 
tems whose Hamiltonian is a sum of kinetic and potential 
energies, the energy surface constraint S(E — iJ(p, q)) 
fixes only the magnitude of the wavevector. The eigen- 
functions therefore appear locally as a sum of randomly 
phased plane waves pointing in arbitrary directions with 
fixed wavevector k. The central limit theorem asserts 
such waves are Gaussian random. An example is shown 
in Fig. ([}]) for a two-degree-of-freedom system where the 
spatial correlations fall off as a Bessel function, Jo(kr). 

If the eigenstates truly possessed these characteristics, 
then a perturbation of the Hamiltonian would have ma- 
trix elements that behaved as Gaussian random vari- 
ables whose variance depended only on the energy sep- 
aration of the two eigenstates, i.e. an energy-ordered, 
banded random matrix. The energy ordering separates 
the weakly interacting states, and therefore only the local 



structure is of importance here. The range of the aver- 
aging carried out in the correlation function is taken to 
be much less than the bandwidth of such a random ma- 
trix. The ultimate statistical expression of this structure 
is embodied in one of the standard Gaussian ensembles 
(GE). We construct a parametrically varying ensemble 

{H(X}} as 

H(X) =H + XHi (4) 

where Hq and Hi are independently chosen GE matrices. 
Note that the sum of two GE matrices is also a GE ma- 
trix which thus satisfies our desire to consider stationary 
statistical properties as A varies. 



FIG. 3. Realization of random wavefunctions in two de- 
grees of freedom. A superposition of 30 plane waves with 
random direction and phase shift, but fixed magnitude of the 
wavevector is shown. 

It is unnecessary to specify the abstract vector space 
of {H(X)} (only the dimensionality of the space) in the 
definition of the ensemble. However, \a) has to be over- 
lapped with the eigenstates, and thus a localized wave 
packet seemingly must be specified. In fact, the specific 
choice is completely irrelevant because the GEs are in- 
variant under the set of transformations that diagonalize 
them. I a) can be taken as any fixed vector in the space 
by invariance. The overlaps and level velocities turn out 
to be independent over the ensemble since diagonalizing 
{Hq} leaves {Hi} invariant and the level velocities are 
equal to the diagonal matrix elements of H-y. With the 
overbar denoting ensemble averaging, 

(p^Mm = ^e(^§t)e =q (5) 

(Ta<TE OaOE 

In fact, it is essential to keep in mind that every choice 
of I a) gives zero correlations within the random matrix 




4 



framework. The existence of even a single \a) in a par- 
ticular system that leads to nonzero correlations violates 
ergodicity. 

It is straightforward to go further and consider the 
mean square fluctuations of C Q (A), 
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where iV is the effective number of states used in the 
energy averaging. Again the level velocities are indepen- 
dent of the eigenvector components. The dEj(X) / dX = 
(J\H±\j) and thus the i ^ j terms vanish due to the in- 
dependence of the diagonal elements of the perturbation 
leaving only the diagonal terms that involve the quan- 
tities that respectively enter the variance of the eigen- 
vector components and the mean square level velocity. 
The final result reflects the equivalence of ensemble and 
spectral averaging in the large- TV limit. Therefore, in er- 
godically behaving systems, C Q (A) = 0± iV -1 / 2 for every 
choice of \a). Fig. (|2|) was made using the orthogonal 
GE. It illustrates a manifestation of ergodicity, i.e. uni- 
versal response of the quantum levels with respect to A 
and democratic behavior of the overlap intensities. 



IV. SEMICLASSICAL DYNAMICS 

We develop a theory based upon semiclassical dynam- 
ics which explains how nonzero overlap correlation co- 
efficients arise out of the localization properties of the 
system. The theory simply reflects the quantum mani- 
festations of finite time correlations in the classical dy- 
namics. In a chaotic system, the classical propagation of 
Pa(p, q) will relax to an ergodic long time average. How- 
ever, wave packet revivals in the corresponding quantum 
system earlier than this relaxation time can occur 
In Heller's original treatment of scars Jig] , he uses argu- 
ments based upon these recurrences which occur at finite 
times to infer localization in the eigenstates. 

In the correlation function, the intensities, p an , weight 
most heavily the level motions of the group of eigenstates 
localized near p a (p, q), if indeed such eigenstates exist. 
If we construct the Hamiltonian as in Eq. (|J) where Hq 
is the unperturbed part, then by first-order perturbation 
theory, the level velocities are the diagonal matrix ele- 
ments of Hi just as in random matrix theory. We showed 
in the previous section that in random matrix theory 



these elements weighted with the intensities are zero cen- 
tered. For a general quantum system the equivalent ex- 
pectation would be fluctuations about the correspond- 
ing classical average of the perturbation over the micro- 
canonical energy surface, S(E — H(p,q)). In this case, 
C a (X) w for all \a). On the other hand, the quantum 
system will fluctuate differently if there is localization in 
the eigenstates. Note that this means some choices of 
| a) will still lead to zero correlations. It only takes one 
statistically significant nonzero result to demonstrate lo- 
calization conclusively, but to obtain a complete picture, 
it is necessary to consider many \a) covering the full en- 
ergy surface. 

We begin by examining the individual components of 
the overlap correlation coefficient, the level velocities and 
intensities. Their ^-dependences are derived and also 
they are shown to be consistent with random matrix the- 
ory as h — > 0. Finally, the weighted level velocities are 
discussed. We give an estimate based upon a semiclas- 
sical theory involving homoclinic orbits for the slope of 
the large intensities. 



A. Level velocities 

In random matrix theory (RMT) level velocities are 
Gaussian distributed as would also be expected of a 
highly chaotic system in the small fi limit. Thus, the 
mean and variance, crfj, give a complete statistical de- 
scription in the limiting case and are the most important 
quantities more generally. Since the purpose of this sec- 
tion is to derive their scaling properties, it is better to 
work with dimensionless quantities. Thus, the dimen- 
sionless variance is defined as a\ = d (E, X)a% where 
d(E, A) is the mean level density which is the reciprocal 
of the mean level spacing. 

We begin by following arguments originally employed 
by Berry and Keating |35| in which they investigated the 
level velocities normalized by the mean level spacing for 
classically chaotic systems with the topology of a ring 
threaded by quantum flux. In order to make the discus- 
sion self contained we will summarize their basic ideas 
using their notation and then extend their results to in- 
clude level velocities for any classically chaotic system. 
More recently, Leboeuf and Sieber Q studied the non- 
universal scaling of the level velocities using a similar 
semiclassical theory. The ^-dependence of the average 
and root mean square level velocities for an arbitrary 
parameter change is derived and is consistent with the 
previous works. 

The smoothed spectral staircase is 



N e {E, A) = ®<i{E — E n (X)) 



(7) 



and taking the derivative with respect to the parameter, 
we obtain 
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dN e (E,X) 
dX 



Y,S e (E - E n (\))?^ (8) 



The quantity e is an energy smoothing term which will 
be taken smaller than the mean level spacing. Our cal- 
culations will use Lorentzian smoothing where 

The energy averaging of Eq. (0) yields 



dN ^\ =- d (E,X)/ dE ^ 



dX I, 



dx 



(10) 



Thus, in order to obtain information about the level ve- 
locities, we will evaluate the spectral staircase. 

The semiclassical construction of the spectral staircase 
is broken into an average part and an oscillating part 



N t (E, X) = N{E, X)+J2 B p{ E i A ) ex P 



S P (E, X) 



x exp 



-eT p (E, X) 



(11) 



The average staircase N(E, X) is the Weyl term and to 
leading order in H is given by 



N(E,X) = ± S 



6(E-H(p,q;X))dpdq (12) 



This simply states that each energy level occupies a vol- 
ume h d in phase space. A change in the phase space 
volume will produce level velocities due to the rescaling. 
We wish to study level velocities created by a change in 
the dynamics, not the rescaling. Hence, without loss of 
generality we will require the phase space volume to re- 
main unchanged, so dN(E, X)/dX = 0. The oscillating 
part of the spectral staircase is a sum over periodic or- 
bits. In general, a perturbation will alter the value of the 
classical actions, S p , the periods, T p , and the amplitudes, 



B p = 



exp(iu p 



2 7 r v /det(M p - 1) 



(13) 



where M p is the stability matrix and v p is the Maslov 
phase index. The summation is most sensitive to the 
changing actions and periods because of the associated 
rapidly oscillating phases, i.e. the division by fi in the 
exponential. Since the energy smoothing term, e, is taken 
smaller than a mean level spacing, it scales at least by 
h d and the derivatives of the period vanish as Ti —* 0. 
Thus, only the derivatives of the actions are considered, 
and the oscillating part of the staircase yields 



dN osc {E,X) 
dX 



\ p 

x exp \i 
x exp 



i dS p {E,X) 

n ox 

S P (E,X) 



h 

eT p (E,X) 



(14) 



It has been shown 
a periodic orbit is 

dS p = 
dX 



that the change in the action for 



r » c?ff(p,q;A) 
dX 



dt 



(15) 



The above integral is over the path of the unperturbed 
orbit and the Hamiltonian can have the form of Eq. (Q) 
where Hq is the unperturbed part. Eq. (|l4| ) can be solved 
without the explicit knowledge of the periodic orbits in 



jlaced by its 
the collec- 



the k — > limit. The quantity dS p /dX is re: 
average. By the principle of uniformity 
tion of every periodic orbit covers all of phase space with 
a uniform distribution. Thus, the time integral can be 
replaced by an integral over phase space upon taking the 
average, 



lim — 



1 /dS* 



dX 



1 



cW(p,q;A) 



V J dX 

xS{E - H(p,q;X))dpdq (16) 



where V is the phase space volume of the energy surface. 
The above treatment of the average is only valid for the 
long orbits, but we may ignore the finite set of short 
orbits in the sum for small enough fi. dH(p, q; X)/dX 
is the perturbation of the system that distorts the en- 
ergy surface. Since the phase space is assumed to remain 
constant, then the average change in the actions of the 
periodic orbits is zero in the limit of summing over all the 
orbits. If only a finite number of orbits are considered, 
corresponding to a finite h, then there might be some 
residual effect of the oscillating part which will cause a 
deviation from RMT. 

Continuing to follow Berry and Keating, the mean 
square of the counting function derivatives can be ex- 
pressed in terms of the level velocities 



fdN t 



\ dX 



E \ n m 

xS e (E-E n (X)) 
xS e {E-E m (X))) E (17) 



For a non-degenerate spectrum, the summation is non- 
zero only if n = m because of the product of the two 
delta functions. Since Lorentzian smoothing is applied, 
then 



for e <C d ■ Thus we have 



'dN e 
dX 



(E,X) 





2ne \\ dX 



(18) 



(19) 



The final result will be independent of e and the type 
of smoothing, i.e. Lorentzian or Gaussian. Using the A 



G 



derivative of Eq. (|Tl|), the dimensionless level velocities 



are 




p p' 



exp <; [T p 



(20) 



The diagonal and off-diagonal contributions are sepa- 
rated, so 



5% 



-2 

a E,di, 



+ a E,off 



(21) 



As 7i — > 0, the phase of the exponential oscillates rapidly 
and averages out to be zero unless S p — Sp>. We will 
assume that this occurs rarely except when p = p' . The 
product (dSp/dX)(dS p > /d\) can take on both positive 
and negative values. This also helps to reduce the con- 
tributions of the off-diagonal terms. For a more complete 
discussion of the diagonal vs. off-diagonal terms see Q. 
We will only present the results for the diagonal terms, 
since the correlations between the actions of different or- 
bits is not known but should not alter the leading h- 
dependence. 

The diagonal contribution is 



This property is independent of the placement of the two 
points, i.e. the two points can lie on the same orbit as long 
as the time between the points increases to infinity. Thus, 
by the central limit theorem, dS p /dX will be Gaussian 
distributed for the sufficiently long periodic orbits. The 
time dependence of Eq. (|23| ) is approximated by a method 
discussed by Bohigas et al. pOfl. They define 



K(E) 



, ^aff(p(0).q(n):Ai ( i/./(p!/).q(/):,\i 



dX 



ox 



dt 

v 

(25) 



which can be evaluated in terms of properties of the per- 
turbation. The variance of the actions in the limit of long 
periods becomes 



(26) 



Applying the Hannay and Ozorio de Almeida sum 
rule (B8[ , the following substitution is made 



(27) 



-2 

a E,diag 




exp 




(22) 

The factor g depends on the symmetries of the system. 
For systems with time-reversal invariance 3 = 2 and with- 
out time-reversal symmetry g = 1. The precise values of 
dSp/dX are specific to each periodic orbit rendering the 
sum difficult to evaluate precisely. A statistical approach 
is possible though which generates a relationship between 
the sum and certain correlation decays. Hence, the quan- 
tity (dSp/dX) 2 in Eq. is replaced by its average, 



T I dH(p(t), q(£); A) 
o \ dX 

8HW)MtJ,X)\ dt , dt 

2 j T / T ( ag(P(t) ' qW;A) 
dH{p(t' + t),(i(t' + t);X) 
OX 




Long orbits increasingly explore the available phase space 
on an ever finer scale. As the time between two points 
in a chaotic system goes to infinity, then they become 
uncorrelated from each other. This is a consequence of 
the mixing property, 



Hence, the diagonal contribution is 



a E,diag 



edg_ 
irh 2 



(2K(E)T) exp 



-2eT 



Jeff 



gK{E)d 



oc K 



-(<*+!) 



(28) 



The variance of the level velocities on the scale of a mean 
spacing grows hT 1 faster than the density of states as the 
semiclassical limit (h — > 0) is approached; see numerical 
tests performed on the stadium in the next section. 

The exact level velocities are perturbation dependent 
and cannot be determined without specific knowledge of 
the system (i.e. the evaluation of K(E)). K(E) is a clas- 
sical quantity that contains dynamical information about 
the periodic orbits, ft should scale as the reciprocal of 
the Lyapunov exponent pl| . Leboeuf and Sieber derived 
K(E) for billiards where the perturbation is a moving 
boundary. In this case K(E) depends upon the autocor- 
relation function and the fluctuations of the number of 
bounces. For maps K(E) is an action velocity diffusion 
{dSp/dX} being Gaussian distributed is 
linked to the level velocities being Gaussian distributed 
as in RMT. If the {dS p /dX} are not Gaussian distributed 
by the Heisenberg time, then one should not expect the 
level velocities to be consistent with RMT; again see the 
stadium results ahead. 



dt'dt (23) coefficient (42 



(/(0)/(% 







(24) 
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B. Overlap intensities 

Now we investigate the overlap intensities and derive 
a semiclassical expression for the Ti-scalmg of the root 
mean square. Eckhardt et al. |Q developed a semi- 
classical theory based on periodic orbits to obtain the 
matrix elements of a sufficiently smooth operator. How- 
ever, the projection operator of interest here, |a)(a|, is 
not smooth on the scale of % for Gaussian wave pack- 
ets. Thus, their stationary phase approximations do not 
apply, in principle, to the oscillating part of the strength 
function. In Berry's work on scars |4J] , he used Gaussian 
smoothing of the Wigner transform of the eigenstates to 
obtain a semiclassical expression for the strength of the 
scars. His approach led to a sum over periodic orbits. We 
will use the energy Green's function similar to Tomsovic 
and Heller in where they derived the autocorrelation 
function using the time Green's function and gave results 
for the strength functions as well. This technique results 
in a connection between the overlap intensities and the 
return dynamics, namely the homoclinic orbits. 

For completeness, we present the smooth part of the 
strength function which is easily obtained from the zero- 
length trajectories, 

S a (E, A) = ^ J A(q, p)6(E - H (q, p))dqdp (29) 

yl(q, p) is the Wigner transform of the Gaussian wave 
packet and is given by 

A(q, p) = 2 d cxp{-(p - p Q ) V/ft 2 - (q - q Q ) 2 /a 2 } 

(30) 

The above results were previously used by Heller [ [f5| in 
the derivation the envelope of the strength function and 
does not contain any information about the dynamics of 
the system. 

The oscillating part of the strength function, on the 
other hand, includes dynamical information, 

S a , osc (E,\) = -^Im J {a\q)G(q,q';E)(q'\a)dqdq' 

(31) 



where 



G(q,q';E) = - 



I -D.. 1 1/2 



xe 



i[Sj(q,q';£;)/a-^7r/2] 



(32) 



is the semiclassical energy Green's function. The above 
sum is over all paths that connect q to q' on a given 
energy surface E. The action is quadratically expanded 
about each reference trajectory; see Appendix ^ for de- 
tails. The initial and final points (q^ and q/) of the 
reference trajectories are obtained by considering the 



evolution of the wave packet. Nearby points will be- 
have similarly for short times. Thus, the phase space 
can be partitioned into connecting areas. As the time 
is increased the number of partitions grow and the size 
of their area shrinks. The reference trajectories are the 
paths that connect the partitions. The autocorrelation 
function in pH has the same form as Appendix ^ where 
the paths that contribute to the saddle points are the 
orbits homoclinic to the centroid of the Gaussian wave 
packet so that q^ and q/ lie on the intersections of the sta- 
ble and unstable manifolds. The result from Appendix |A| 
is 



S a ,osc{E) - ^^Re^ 



det A 21 
det A 



1/2 



(33) 



1/2 



\q [N) \\q' [N) \ 



where 



1 i 
/i(q/,qi) = exp \ -b • A" 1 • b - -p a ■ (q f - q*) (34) 



(q/ - qq) 2 (q» - qq) 2 



2CT 2 



2ct 2 



The function /j(q/, q%) in the above equation is a damp- 
ing term which depends on the end points of the ho- 
moclinic orbits. Only orbits which approach the center 
of the Gaussian wave packet in phase space will con- 
tribute to the sum. The time derivatives of the paral- 
lel coordinates are evaluated at the saddle points which 
are near the centroid of the Gaussian, so we may set 
« Ig'^l ss |p Q |/?7i. The sum over homoclinic 
orbits used for the autocorrelation function in |2l[] con- 
verged well to the discrete quantum strength function 
when only those orbits whose period did not exceed the 
Heisenberg time, (th = 2tt Tid(E, A)), were included. As 
happened with the periodic orbits and the level veloci- 
ties, in order to evaluate Eq. (33) the homoclinic orbits 
and their stabilities must be computed rendering the sum 
tedious to evaluate precisely as done in 

By taking a statistical approach we can gain some in- 
sight into the workings of this summation. The variance 
of the intensities are obtained by a similar fashion as the 
level velocities. Using Eq. (|j) and Eq. (|lj), we have 



<^, OS c(^A)>. 



2ttc 



(35) 



Since the square of the strength function is a product 
of two delta functions, an energy smoothing term is re- 
quired. After making the diagonal approximation, we 
obtain 



a,diag 




(36) 
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A classical sum rule is applied to the above sum for 
special cases including two-dimensional systems; see Ap- 
pendix H for the details. Thus, 



2egm 2 a'- 
fi?d\p a \' 

2 2 

gm z a z 
M\p a \ 2 



J exp{-2eT/h}dT 



(37) 



1 II 

Setting a oc h ' which shrinks the momentum and posi- 
tion uncertainties similarly, the ?i-scaling of a 2 a diag is Ti d \ 
see numerical tests of the stadium in the next section. 

Assuming that the amplitudes of the wavefunctions 
are Gaussian random, then the RMT result for strength 
functions is a Porter-Thomas distribution which has a 
variance that is proportional to the square of its aver- 
age. The average strength function, Eq. (p9|), scales as 
{<j/Ti) d for Hamiltonians which can be locally expanded 

1/2 

as a quadratic. Therefore, with cr oc h 1 the variance of 
the strength function, Eq. (|37|), scales as the square of 
the average and is consistent with RMT. 



C. Weighted level velocities 

A semiclassical treatment of the overlap correlation co- 
efficient defined in Eq. (|J) is now developed. As stated 
in the introduction, the companion paper p3| presents 
the semiclassical theory for maps. We stress that in the 
preceding subsections and in what follows is for conser- 
vative Hamiltonian systems. Here, the ^-dependence of 
the average overlap correlation coefficient is established 
and a semiclassical argument for the existence of nonzero 
correlations is presented. 



1. Actions of homoclinic orbits 

To calculate the overlap correlation coefficient, the rate 
of change of the actions for homoclinic orbits will be nec- 
essary. As discussed earlier, this was accomplished for 
periodic orbits |37). We extend these results to include 
the actions of homoclinic orbits. Homoclinic orbits have 
infinite periods causing their actions to become infinite. 
We are interested in the limiting difference of the action, 
<S,- , between the j th homoclinic orbit and repetitions of 
its corresponding periodic orbit p. The difference is finite 
and is equal to the area bounded by the stable and un- 
stable manifolds with intersection at the j homoclinic 
point in a Poincare map. Sj provides information about 
the additional phase gathered by the homoclinic orbit. 
The action of the j th homoclinic orbit as n — > oo in the 
time interval (—nT p ,nT p ) is 



:(p) 



2nS v + S 



(p) 



(38) 



where T p and S p are the period and action of the periodic 
orbit, respectively. As a consequence of the Birkhoff- 
Moser theorem , if the Poincare map is invertible and 
analytic, then there exist infinite families of periodic or- 
bits that accumulate on a homoclinic orbit. It is thus 
possible to estimate the action of the homoclinic orbit 
by these periodic orbits whose action is given by |47| 



= nS p + Sf J -s 



(p) » 



(39) 



(v) 

where s^ j is the difference in action between a path de- 
fined by Sj P ^ along the stable and unstable manifolds and 

the path of the new periodic orbit in a Poincare map. 

(p) (p) 
s„ ij depends exponentially on n, so as n — > oo, a^nj 

approaches the action of the homoclinic orbit. Thus, in 



the limit of large n, 5- is approximated by the differ- 



ence between two periodic orbits (i.e. S 



(v) 



Hence, the change in S 
calculated as in 11371, 



(p) 



Ap) 

n,3 



nS v 



due to a small perturbation is 



AS 



(p) 



-AA 



-nAA 



dH(p,q; A) 

bo d\ 

dg(p,q;A) 
dX 



dt 



dt + 0(X 2 



(40) 



where the integrals are over the unperturbed periodic or- 
bits. The differences, s~j, can be made smaller than the 
second order term in Eq. ( [Io"| ) by taking n large enough. 
Interchanging the order of integration and differentiation, 
the integrals reduce to the unperturbed energy times the 
derivative of the orbit period with respect to the param- 
eter, 



AS. 



(p) 



-AXE, 



d_ 

p dx 



- nTr 



The orbit period T can be expressed as dS/dE. 
the difference of the two periods as n — > oo is 



T 



da 



(p) 

n ,3 



8E 



dS n dS. 



dE 



(p) 



dE 



(41) 



Thus, 



(42) 



Hence, 



AS 



(p) 



-AA- 



dX 



' dS^ 

E, 3 



-AA 



p dE J 

dH(p,q;X) 
dX 



dt 



(43) 



Note that the integral is over the unperturbed path along 
the stable and unstable manifolds. For zero correlations, 
AiS-^ must be "randomly" distributed about zero. 

If enough time is allowed, then for ergodic systems 
the set of all homoclinic orbits for a given energy will 
come arbitrarily close to any point in phase space on 
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that energy surface. Thus, an integral over phase space 
on the original energy surface can be substituted for the 
time integral, 



1 



7<P)\ 



lim - ( AS) > 



J S(E-H(p,diO)) 

dpdq (44) 



-AA 

,dff(p,q;A) 



dX 

Since the density of states are kept constant, the pertur- 
bations fluctuate about zero and the integral vanishes. 
Hence, the average change in the actions will be zero. 
The mean square fluctuations of the actions, 



(ASf) ) =(AA) 

3 



JO 

dH(p',q!;\) 
dX 



dH (p,q;A) 
dX 

dtdt' (45) 



approaches a Gaussian distribution by the central limit 
theorem via the same reasoning as that for periodic or- 
bits. Again we can define Kh om {E) as in Eq. (25|), except 
now the average is over homoclinic orbits and instead of 
integrating to infinity we only integrate to the Heisen- 
berg time to be consistent with the range of the sum in 
Eq. ([33l). It is the short time dynamics that dominate. 
Long time correlations will average to zero by the mixing 
property (Eq. (E3)). Thus, the variance of the actions 
becomes 




2K hom (E)T 



(46) 



Khom(E) will approach K{E) in the semiclassical limit 
(t h -> 00). 



Lorentzian smoothing was again employed, Eq. (|18|), and 
we've defined C a (X) to be the numerator of the overlap 
correlation coefficient (without the division of the rms 
level velocities and intensities) . By the definition of the 
overlap correlation coefficient only the oscillating part 
of the level velocities and the intensities are considered. 
Using the derivative with respect to lambda of Eq. (|T] 
and Eq. (f^) the numerator becomes 



- 2ire / v^v^ mcrBr, ( detA 21 

C a (X) = — (Rc^^ 

j v 



1/2 



d \ itTi 
\ 3 p 



detA 



(48) 



\p a \ 



dX ) A 



Because of the rapidly oscillating phases, the energy aver- 
aging will result in zero unless Sj ~ S p . As stated earlier, 
for every homoclinic orbit there is a periodic orbit that 
comes infinitesimally close to it. The same periodic or- 
bit's action can be nearly equal to the actions of different 
segments of the same homoclinic orbit. Thus, a diagonal 
approximation is used for the homoclinic segments 



C„(A) 



2-Ktg j n ^ maBj I det A 21 



1/2 



^ ReV- r 
d \ 4^ irh 2 



9S A /j(q/.qi) 



-2eTj/h 



detA 



dX 



[Po 



(49) 



Upon applying the sum rule f or t wo-dimensional systems 
and other special cases, Eq. (B8), we have 



4(A) 



egma 1 



7td\p 



h 2 jo 



3 -2eT/ft 



(50) 



2. Overlap intensity-level velocity correlation coefficient 

In the previous two subsections, we have examined the 
pieces that constitute the overlap correlation coefficient. 
The semiclassical theories of the level velocities and the 
intensities are now combined to construct a semiclassical 
theory for the weighted level velocities. The numerator 
of the overlap correlation coefficient is proportional to 
the energy averaged product of the intensities and level 
velocities, 



S a (E,X) 



dN{E,X) 
dX 



= (^2^2Pan(X) 



dE m 
dX 



xS e (E-E n )5 e (E-E m )) E 
d / dE n 



(47) 



The changes in action of the homoclinic excursions are 
now weighted by the /j(q/, qi)'s. Without the additional 
weighting the average in the changes in the action would 
be zero for all positions of the Gaussian wave packet. 

In 0|, a heuristic argument for the direction of the 
weighted level velocities was given. The argument ba- 
sically states that the energy surface changes with the 
parameter such that the action changes are minimized. 
Eq. j5C| ) differs from [B2| in that the proper weightings, 
fjififiQi)-) of the homoclinic orbits are derived here, and 
the action changes are not correlated with the inverse 
periods. Also, in the homoclinic orbits were strictly 
cut-off at the Heisenberg time whereas here there is an 
exponential decay on the order of the Heisenberg time 
with the energy smoothing term, e, equal to Ti/th [ [35"[ . 
One reason that j22] reported such good results is that 
since the number of homoclinic segments proliferate ex- 
ponentially, most of the included segments occurred near 
the Heisenberg time and the expression in Eq. ( poj ) is 
divided by the Heisenberg time (i. e. multiplied by e). 
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As Ti — * (tjj — > oo), the integral in Eq. (50) would 
be dominated by the subset of {dSj/dX} associated with 
very long orbits and would decouple from the weightings. 
For small enough Ti, as previously stated, the dSj/dX 
for these orbits would approach a zero-centered Gaussian 
density, and the integral would vanish. In other words, 
we could use arguments analogous to those underlying 
Eq. (O) to write 



C a (X) 



-gma 



7rd\p a \ 



dpdq ag( ^ q;A) J(£;-g(p,q;A)) (51) 



where F is the phase space average of the weightings. 
Note that the phase space average of <9if (p, q; A) / d\ van- 
ishes excluding an irrelevant drift of levels, so the RMT 
prediction of C a (X) is recovered for h small enough. 

The leading order in h correction to this is more dif- 
ficult to ascertain. H enters into the exponential in the 
integral for the energy smoothing, but not for the clas- 
sical decay of the action changes. Upon taking the in- 
tegral, this yields two competing terms for the 7i-scaling 
which may depend upon the region of phase space the 
correlation is taken in. The numerics also show a large 
fluctuation of the scaling in the stadium (see the next 
section). 



V. STADIUM BILLIARD 



altering the radii of the endcaps, R, to keep the area 
of the stadium a constant. Throughout this section, the 
level velocities and intensities are evaluated for a stadium 
with A = R = 1. For billiards the average number of 
states below a given energy, E, is approximately N(E) « 
mAE /2irfi 2 where A is the area of the billiard. This is 
the first term in an asymptotic series in powers of h. 
The density of the states, dN/dE, is then a constant 
not depending on A to the lowest power of h if the area 
remains the same. 
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FIG. 4. Birkhoff coordinates for the stadium billiard. The 
position coordinate can be taken to start anywhere along the 
perimeter. Here we have chosen the origin to be the middle 
of the right semicircle. 

In this section the semiclassical theories just presented 
and the numerical results from the stadium billiard are 
compared. The stadium billiard, which was proven by 
Bunimovich E§] to be classically chaotic, has become a 
paradigm for studies of quantum chaos. It is defined as a 
two-dimensional infinite well with the shape pictured in 
Fig. (|j). We continuously vary the side length, 2 A, while 



FIG. 5. Distribution of the level velocities for a stadium 
billiard. The solid line is the lowest energy range, the dotted 
line is the middle energy range and the dash-dot line is the 
highest energy range. The RMT result is given by the dashed 
line. The level velocities have been rescaled to zero mean and 
unit variance. 

We will examine three different energy regimes for the 
stadium. Since billiards are scaling systems, this will 
correspond to three different values of h. The energy 
regimes are separated by a factor of four in energy or 
conversely a factor of one half in h. Twice as many states 
are taken in each successive energy regime so that the 
averages will incorporate the same relative size interval 
in energy as h is decreased. This corresponds to the 
increase in the density of states for varying Ti. 

The distributions of the level velocities for all three 
energy regimes are shown in Fig. (|^) along with the ran- 
dom matrix theory prediction. The skewness occurs be- 
cause of a class of marginally stable orbits in the sta- 
dium. These orbits are the bouncing ball orbits which 
only strike the straight edges. Their contribution do not 
seem to decrease as the semiclassical limit is approached 
though they should once h is sufficiently small. There 
is no clear trend for the level velocity distribution to 
approach Gaussian behavior. The root mean square of 
the level velocities also deviates from our calculations 
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of the ^-scaling in Section IV (Fig. (g)). This is again 
explained by the bouncing ball orbits whose effects are 
missing from the trace formula. Quantizing only these 
orbits using WKB yields a dimensionless level velocity 
scaling of K~ , while the trace formula gives a scaling of 
7j-3/2 The numerical results give a scaling of approx- 
imately ?i -1 ' 8 which lies in between the two suggesting 
that the marginally stable orbits significantly effect the 
level velocities. 
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FIG. 6. Root mean square of the level velocities as a func- 
tion of 1/h. The solid line is the theoretical value from Sec- 
tion 



IV 



and the dashed line is the WKB results for the bounc- 
ing ball motion. The best fit line through the stadium results 
is the dotted line. 

To study the intensities, the eigenstates must first be 
constructed. Bogomolny's transfer operator method |49[] 
was used to find the eigenstates. This method uses a 
(d — 1) dimensional surface of section. A convenient 
choice is the boundary of the stadium (Fig. (0)). The 
generation of a full phase space picture of the stadium 
would otherwise require four dimensions, two positions 
and two momenta. The position coordinate is measured 
along the perimeter and the momentum coordinate is de- 
fined by cos 9. The classical dynamics have a quantum 
analog that uses source points on the boundary. Thus, 
all of the eigenfunction's localization behavior can be ex- 
plored using wave packets defined in these coordinates. 
A coherent state on the boundary is a one-dimensional 
Gaussian wave packet; see the lower figure in Fig. (Q). 
The corresponding wave packet in the interior of the 
stadium can be generated by a Green's function and is 
shown in the upper figure of Fig. (0). For billiards the 
Green's function is proportional to a zero" 1 order Hankcl 
function of the first kind, H^\kr) /2i% 2 . The centroid of 
the Gaussian wave packet is moved along the boundary 



and its momentum is changed according to the Birkhoff 
coordinate system. Thus, the entire phase space of the 
stadium is explored. 

The results for the average and the standard deviation 
of the intensities using Birkhoff coordinates are shown in 
Figs. (||) and (j^), respectively. The average is flat ex- 
cept for peaks associated with the two symmetry lines 
of the stadium. The eigenstates used here were even- 
even states, so there is twice the intensity along the two 
symmetry lines that bisect the end caps and the straight 
edges. The standard deviation has two large peaks cen- 
tered around the bouncing ball orbits. The rest of the 
figure is relatively flat with a few small bumps. Random 
matrix theory would predict this to be a flat figure with 
small oscillations. The marginal stability of the bouncing 
ball orbits can be seen but no other feature of the sta- 
dium, except for the horizontal bounce, is picked out by 
looking at the intensities. Fig. ([To]) shows the ?i-scaling of 
the root mean square for the intensities where the wave 
packet is placed on various periodic orbits. The theory 
from Section IV predicts a smaller scaling than the nu- 



merical results of the stadium. 




Q 



FIG. 7. The lower figure is Gaussian wave packet on the 
boundary. The upper figure corresponds to the wave packet 
in the interior of the stadium. 
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FIG. 8. Average overlap intensity for a Gaussian wave 
packet denned in the Birkhoff coordinates of the stadium. 



ists phase space localization for most of the states. The 
large positive values of the correlation coefficient in the 
center of the figure again correspond to the bouncing 
ball states. Classically, this area of phase space is dif- 
ficult to enter and leave. Hence, the localization is ex- 
pected to be stronger for this area of phase space. The 
area beneath the peaks is several standard deviations 
(N^ 1 / 2 = (114) -1 / 2 w 0.09) away from zero as predicted 
by random matrix theory. Thus, phase space localiza- 
tion is also occurring in this region. The point exactly 
in between the peaks is the point in phase space asso- 
ciated with the horizontal bounce. The series of smaller 
peaks leading up to the large peaks are the gateways into 
the vertical bouncing ball area. Fig. (|l^) is a plot of the 
orbits corresponding to these peaks. They are periodic 
orbits which only strike the endcaps twice and become 
almost vertical. Orbits must pass through these regions 
in order to enter or exit the vertical bounce states. 



The heights of the bouncing ball peaks can be approx- 
imated by quantizing the rectangular region of the sta- 
dium. The intensities obtained from this calculation are 
weighted by the ratio of the density of the bouncing ball 
states to the total density of states. The Gaussian 
wave packet is placed in the center of the straight edge 
and the middle energy regime is used. The results of this 
approximation are 80.6 and 320.5 for the average inten- 
sity and rms intensity, respectively, compared to 80.7 and 
386.3 for the numerical calculations of the stadium. 




FIG. 9. Root mean square of the overlap intensity for a 
Gaussian wave packet defined in the Birkhoff coordinates of 
the stadium. 
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FIG. 10. Root mean square of the overlap intensities as a 
function of 1/h. The wave packet was placed on the horizontal 
(solid circles), V (solid squares), diamond (solid triangles), 
rectangle (open circles) and bow tie (open squares) orbits. 
The solid line is the theorectial value of h 1 from Section [iv| 
and the dotted line is the best fit of the stadium results which 
is Ti 1 ' 3 . The intensities have been rescaled at Ti = 1 so that 
they occupy the same area of the plot. 



Random matrix theory suggests that the correlation 
coefficient for a generic chaotic system should result in 
zero. On the other hand, using the correlation coefficient 
for the stadium in Birkhoff coordinates, we found that 
some of the states gave nonzero correlations, Fig. (11). 
In fact, large correlations are found for nearly all the 
states in the stadium billiard which means that there ex- 
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FIG. 11. Overlap correlation coefficient for the stadium us- 
ing Birkhoff coordinates. The energy range of the averaging 
is 2200-2600 where h = m = 1. 

As the energy of the system is increased (i. e. Ti is de- 
creased), the results of the correlation function remain 
qualitatively the same, Fig. (|l3|). All the peaks and val- 
leys stay in the same place. The numerical results of the 
overlap correlation coefficient fluctuate depending upon 
the area of phase space being considered. This is consis- 
tent with the the semiclassical theory in Section IV. More 



details of the system are explored as fi is decreased, since 
the phase space is divided into finer areas. Thus, more 
detailed information about the phase space localization 
of the system is observed in the overlap correlation coef- 
ficient at smaller values of Ti. 




Q 



FIG. 12. Trajectories corresponding to the peaks leading 
up to the bouncing ball orbits. The lower figure is a contour 
plot of Fig. The solid circles correspond to the bounce 

points of the trajectories. Geometric and time-reversal sym- 
metries were also included. 



VI. CONCLUSIONS 

We have shown that intensity weighted level veloci- 
ties are a good measure of the localization properties for 
chaotic systems. They are far more sensitive to local- 
ization than similarly weighted level curvatures (which 
are closely related to level statistics). Thus, a system 
can be RMT-like, yet the eigenstates are not behaving 
ergodically (as RMT predicts). 

The stadium eigenstates show a great deal of localiza- 
tion. Not only are the vertical bouncing ball orbits pre- 
dicted by the measure, but also other orbits. The over- 
lap correlation coefficient is very parameter dependent. 
Choosing a different parameter to vary would highlight 
other sets of orbits depending on how strong the pertur- 
bation effects those orbits. The degree of localization can 
be predicted by the return dynamics. In a chaotic system, 
all the return dynamics can be organized by the homo- 
clinic orbits. The manner in which a chaotic system's 
eigenstates approach ergodicity as % — > will depend on 
a new time scale, i.e. that required for the homoclinic 
excursions to explore the available phase space fully. 




FIG. 13. The same as Fig. ([n]) except for a higher energy 
range of 9200-10000 where h = m = 1. Note the finer struc- 
ture of the various peaks. 

Parametrically varied data exist that can be analyzed 
in this way. In the Coulomb-blockade conductance data 
to the extent that the resonance energy variations are re- 
lated to a single particle level velocity (minus a constant 
charging energy and absent residual interaction effects) 
should show correlations. We mention also that the mi- 
crowave cavity data can be studied with even more flex- 
ibility since they have measured the eigenstates and can 
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therefore meticulously study a wide range of \a) to get a 
complete picture of the eigenstate localization properties. 

Finally, this analysis could be applied in a very fruitful 
way to near-integrable and mixed phase space systems. 
In these cases, standard random matrix theory would 
not give the zero" 1 order statistical expectation, but the 
localization would still be determined by the return dy- 
namics in the semiclassical approximation. 
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APPENDIX A: GAUSSIAN INTEGRATION 

Inserting Eq. ( |32| ) into Eq. (|3l|), the strength function 
involves two TV-dimensional integrals where N is the sys- 
tem's number of degrees of freedom, 



Sa - osc ^ - n Im in(2nin)v-w (™ 2 



d/2 



x l^2\D s \ 1 ^exp{-ip a -(q-q')/Ti 

j 

-(q-q Q ) 2 /2a 2 -(q'-q a ) 2 /2<7 2 
+iS j (q, q'; E)/h - ivft /2}dqdq' (Al) 

To evaluate the integrals over q and q', the action is 
quadratically expanded about the points q/ and q^, 

Sj(q.,q!;E) = Sj(qf,qi;E)+p f • (q - q/) - p, ■ (q' - q*) 



1 N 

IE 



(q(i)- q f)( q W- q f) 



dp 



(g'W-^J^W-gW) (A2) 



+ 2 



dp 



(i) 



It is useful to define the vector 

z = (zi, ...,z N ,z[,.. .,z' N ) 

where 

Thus, the integrals become 



(A3) 



(A4) 



-Im 



7T in{2mh)( d - 1 )/ 2 



d/2 



X 



jY,\D s \ 1/2 <J 2d 

3 



x exp{— z • A • z — b • z + c}dz (A5) 
where A is composed of four iV-dimensional matrices 



A 11 A 12 
A 21 A 22 



(A6) 



and 



( f N) > 



with 



and 



b = (iSp^p — Sq^p , . . . , iSp^ — Sq 
i6pV-6qW,...,i5pW-6qW) 

Spf=(p$-pf)<r/h 

5qf={q$-qf)/o 
5qf ) =($-qf ) )/a 



1 1 

c = -S j {q fl q l ;E) - -p a ■ (q f - q*) 
(q/ - q Q ) 2 " Qa) 2 iVji 



(A7) 



(A8) 



2a 2 



2a 2 



(A9) 



The matrix A can be expressed in terms of the stab ility 
matrix, M, where M has the same form as Eq. ( [Aq 1 ) 



(A10) 



Thus, 



li 



2 2ft I flgW 



5g,b I 

2 



cr 2 Ef^.cof (m 2 *) 



27) 



det M 21 



2 




A 12 = _ 

ab 2Ti \dq'W 



M 11 (M 21 ) -1 

(a 
/ 



- / dp {a) 



cof 



2h det M 21 



= ^((M^) 

A 21 = 
ab 2h ydqW J 



21^-l^T 



If 
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1(7 



2 Cof 



2ft det M 21 



A 22 



1(7 

2ft 

8g,b 

2 




det M 21 



la 



{M 21 )- 1 M 2 



(All) 



_ I 

~~ 2 ~ 2ft 

where rn^fc are elements of the stability matrix and 
cof (M 2 ^) is the signed minor of M 2 ^. D s is a deter- 
minate involving second derivatives of the actions, 



d-'S 



'q9q' 6ck)E 

A 2 C C 



dEdq' dE 1 



|9 (JV) II<?' (JV) | 
1 

A 21 



-d 2 S 



<9qdq' 



2/7 



(d-i) 



(A12) 



The tildes in the determinants in the above equation are 
used to exclude the Ath coordinate. To obtain the above 
result q( N ^ and q'( N > are chosen to be locally oriented 
along the trajectory where the dots indicate time deriva- 
tives. Since A is a symmetric matrix, the result for a 
general Gaussian integral is used and, hence, the strength 
function becomes 



^ OSc(jB) ~ 7T ^ft^ft)^ 1 )/ 2 



d/2 



2ft 



(d-l)/2 




can be reduced to determinants of (A x A) matri- 
ces by using the relation ((M 21 ) _1 ) T = -M 12 + 
M n (M 21 ) -1 M 22 and some row and column manipula- 
tions Ell], so that 



x exp 



where the time derivatives are evaluated at the saddle 
points. 

APPENDIX B: SUM RULE FOR THE 
STRENGTH FUNCTION 

The determinant of the (2 A x 2 A) matrix A, 

£((m 21 )-T 



I _ ^M U (M 21 ) -1 



det A = det 2 2 % , 



2ft 

i _ ^(M 21 )-!M 22 



(Bl) 



det A = det 



4h 



M 1 



M 



Am 21 

a 1 



ft 



/det(M 21 ) (B2) 



The coordinates parallel to the trajectory do not mix 
with the transverse coordinates, since a point on an or- 
bit will remain on that particular orbit. Thus, the Ath 
rows and columns of the individual matrices in the above 
expression are zero except for the (A, A) elements. 

It is convenient to re-express the submatrices of the 
stability matrix in terms of the Lyapunov exponents. Let 
{Ai} be the set of Lyapunov exponents whose real part 
is positive ordered such that Ai > A2 > ••• > \n-i- 
The Lyapunov exponents along the parallel coordinate 
are zero and we will only work with the reduced (2(A — 
1) x 2(A — 1)) stability matrix in what follows. Let A 
be the diagonal matrix of the eigenvalues of the reduced 
stability matrix, 



,Xit 













-Ait 



V 
















-Ajv-it 



(B3) 



Thus, by a similarity transform the reduced stability ma- 
trix be can written in terms of the Lyapunov exponents, 
i.e. 



M = LAL 



(B4) 



Hence, each of the elements of the stability matrix can 
be written as 



(A13) 



JV 
k 



-An,* 



(B5) 



where a,-? and are linear combinations of the el- 
ements of the L and L _1 matrices. Because in gen- 
eral chaotic systems A » 0, the £>y s may be omitted 
without seriously effecting the above sum. All the de- 
terminants including the numerator and denominator of 
det A as well as det A 21 , thus, will involve products of 
Eq. (B5). The homoclinic orbits in the sum begin and 
end at the intersections of the stable and unstable mani- 
folds near the Gaussian centroid. Since neither manifolds 
may cross themselves, then in the vicinity of the Gaus- 
sian centroid the branches of each manifold are nearly 
parallel to themselves. Thus, to an excellent approxima- 
tion, the same similarity transformation will diagonalizc 



1G 



the stability matrix for each individual orbit, regardless 
of the period. Consequently, the elements of L and L _1 
are period independent. 

Connections can be made between the determinants 
and the Kolmogorov-Sinai entropy. The Kolmogorov- 
Sinai entropy, h^s, of a system can be expressed using 
Pesin's Theorem as the sum of the Lyapunov exponents 
with positive real part, 



N-l 



IKS 



(B6) 



If there is no mixing between the different coordinates, 
then the individual matrices M 11 , M 12 , M 21 and M 22 
are diagonal. Thus, each matrix element depends only 
upon one Lyapunov exponent and the determinants are 
proportional to exp(— h^gt). This is the case for two di- 
mensional systems where the parallel and perpendicular 
coordinates in the stability matrix separate as mentioned 
above. Hence, we have 



det A 21 



det A 



oc exp(-hKst) 



(B7) 



Unlike the periodic orbits, the homoclinic sum is over 
segments of the orbits. The number of homoclinic points 
will proliferate exponentially at the same rate as the 
fixed points in the neighborhood which is proportional 
to exp(/i^T) where hr is the topological entropy. This 
is demonstrated by examining the partitioning of the 
phase space mentioned in Sect. IV B which has expo- 
nential growth. The partitioning reflects the symbolic 
dynamics of the system. The symbolic code uniquely 
describes each orbit so that amount of code (partitions) 
cannot grow faster than the number of periodic points, 
since each code (partition) cannot represent more than 
one periodic point. 

Finally, the sum rule is obtained by setting the topo- 
logical entropy and Kolmogorov-Sinai entropy equal to 
each other. Then, for the special case of no mixing in the 
stability matrix as mentioned above the combination of 
the amplitudes and the number of orbits yields 



E 



det A 21 



det A 



dT- 



(B8) 
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